Experimental observation of quantum entanglement in low dimensional spin systems 
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We report macroscopic magnetic measurements carried out in order to detect and characterize 
field-induced quantum entanglement in low dimensional spin systems. We analyze the pyroborate 
MgMnE^Os and the and the warwickite MgTiOBOa, systems with spin 5/2 and 1/2 respectively. By 
using the magnetic susceptibility as an entanglement witness we are able to quantify entanglement 
as a function of temperature and magnetic field. In addition, we experimentally distinguish for the 
first time a random singlet phase from a Griffiths phase. This analysis opens the possibility of a 
more detailed characterization of low dimensional materials. 
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Since the development of quantum mechanics, entan- 
glement has been a subject of great interest. Lately, this 
is mainly due to the importance of entanglement in quan- 
tum information and computation. As a consequence, a 
great effort has been made to detect and quantify entan- 
glement in quantum systems Q. In addition, quantum 
spin chains, a class of systems well known in solid state 
physics, began to be studied in the framework of quan- 
tum information theory; there are proposals for the use of 
such systems in quantum computation Q . Naturally, en- 
tanglement in interacting spin chains acquired relevance 
in the QI community. Therefore, there has been a special 
effort in understandingand quantifying quantum entan- 
glement in solid-state 0, 0, 0, . At the same time, the 
condensed matter community has begun to notice that 
entanglement may play a crucial role in the properties of 
different materials [?|. 

It is a difficult task to determine experimentally if 
a state is entangled or not. One of the new promis- 
ing methods for entanglement detection is the use of 
an entanglement witnesses (EW). EW are observables 
which have negative expectation values for entangled 
states. Magnetic susceptibility was recently proposed as 
an EW and some old experimental results were re- 
analyzed wthin this new framework 

It has been known for a long time that entanglement 
appears in quantum spin chains, like the spin 1/2 Heisen- 
berg model. The disordered spin 1/2 one-dimensional 
Heisenberg model, for example, presents a random sin- 
glet phase (RSP), where singlets of pairs of arbitrarily 
distant spins are formed @ . Although entanglement was 
already known to exist in such chains, it had not been 
quantified theoretically until this decade ■ A previous 
study of a diluted magnetic material Q has shown the 
importance of entanglement, but to our knowledge, this is 
the first experimental measurement of quantum entangle- 
ment in a magnetic material. As representative systems, 
we analyze the pyroborate MgMnE^Os 0,0] and the 



warwickite MgTiOBOa 0|, two quasi-one dimensional 
disordered spin compounds with previously known mag- 
netic and thermodynamic properties that suggest the ex- 
istence of either a RSP or a Griffiths phase (GP) 0] at 
low temperatures. 

There are no experimental studies on random mag- 
netic chains which discriminate these two phases. In 
this Letter, from a detailed analysis (if magnetic measure- 
ments, we show unambiguously the existence of a RSP 
in MgTiOBOa, which is expected for a spin-1/2 random 
exchange Heisenberg antiferromagnetic chain (REHAC). 
In addition, our study of MgMnB2 05 provides experi- 
mental evidence for the existence of a Griffiths phase in 
a low dimensional system with S > 1/2. 

Addressing the entanglement properties of these ran- 
dom systems, there is also a clear distinction between the 
RSP and the Griffiths phase. For the former, entangle- 
ment is well characterized and has been shown to scale 
with the logarithm of the size 0,0] . For the latter there 
is no theoretical study of how entanglement behaves. 

We make use of magnetization and ac susceptibility 
measurements as a function of temperature and applied 
magnetic field to detect and quantify entanglement by 
using the susceptibility as an entanglement witness |l(. 
First we show that both systems present entanglement at 
low temperatures with no applied field. Next, we analyze 
the ac. susceptibility and magnetization as function of 
field for different temperatures and we quantify the vari- 
ation of entanglement as a function of applied field. We 
observe that entanglement increases for increasing mag- 
netic fields in a region of the B x T diagram. This unusual 
behavior was suggested by Arnesen et al. and called 
magnetic entanglement. 

In both pyroborate MgMnB20s and warwickite 
MgTiOBOs there are ribbons formed by oxygen octa- 
hedra sharing edges. These octahedra give rise to four 
columns, along the ribbons, whose centers define a tri- 
angular lattice and two different crystallographic sites 
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for metals: one in the central columns and another in 
the border ones. In the pyroborate such columns do not 
touch and both metal sites are equally occupied by the 
two metal ions [lij . In the warwickite, on the contrary, 
the columns do touch and the metal sites are probably 
occupied as in MgScOBOa: 76 % of the internal sites oc- 
cupied by the transition metal and 24 % by the alkaline- 
earth metal. The sites on the border columns have the 
opposite occupancy [l7j . 

The pyroborate powder was obtained from grinded sin- 
gle crystals, and the warwickite powder was directly ob- 
tained from its synthesis. The warwickite sample was 
analyzed through X-Ray diffractometry; it has been ver- 
ified that the material was well crystallized and that its 
purity was 97.72 % , as evaluated by the method of Lut- 
terotti et al. 01- The more abundant impurity was the 
non-magnetic MgTiOs. More details on sample prepara- 
tion were previously published |l2lll3j . Dc magnetization 
and ac susceptibility measurements were performed with 
a commercial apparatus (Quantum Design PPMS). 

In MgTiOBC>3, evidence for a RSP-like behaviour was 
previously obtained from specific heat C(T), susceptibil- 
ity x(T) and magnetization m(H) measurements 
These quantities exhibit the characteristic power law be- 
havior associated with a RSP, x(T) oc T~ a , down to the 
lowest measured temperature. In this system the mag- 
netic ion Ti 3+ has spin S=l/2, and due to the negligi- 
ble magnetic anisotropy this material is well described 
by a spin-1/2 (REHAC). The physical behavior is con- 
trolled by an infinite randomness fixed point indepen- 
dent of the amount of disorder. On the other hand in 
the MgMnB205 pyroborate, the magnetic ion Mn 2+ is a 
spin 5/2, S state ion. The phase diagram of a REHAC 
with S > 1/2 is not a trivial one. For weak disorder these 
systems pre sent GP, while only for strong disorder a RSP 
appears 

In figure 2] panels (a) and (b), we show the ac mag- 
netic susceptibility as a function of temperature for 
MgTiOB0 3 and MgMnB 2 5 respectively. Both sys- 
tems have a sub- Curie regime at low temperatures. 
MgMnB205 acquires a Curie-like temperature depen- 
dence around 50 K. On the other hand, MgTiOB03 
presents a sub-Curie susceptibility even at room temper- 
ature. It is known that both systems have a susceptibility 
which behaves as x(T) cx T~ Q , although the temperature 
dependence of a was not further analyzed. 

These two different phases should be distinguished ex- 
perimentally by the temperature dependence of the expo- 
nent a. The GP is characterized by a constant value of a. 
For the RSP, we should expect a low-temperature suscep- 



tibility following x(T) 



Tln a (f2 /T) 



[2JJ, which is equiv- 



alent to a(T) = 1 — ln ( T ° 0o ) ) where a is a constant |2l|. 
So, the RSP is characterized by a slowly varying a(T). 

Following Hirsch [2l| , we analyze the data by redefin- 
ing a — —d(ln(x))/d(\n(T)) and extract the tempera- 
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FIG. 1: Up: Magnetic susceptibility versus temperature, (a) 
Experimental data for MgTiOBOs (open circles) and fitting 
using the susceptibility expression of a RSP (solid line). On 
panel (b) , the experimental data for MgTiOBOs (open circles) 
and for high temperatures a Curie- Weiss fiiting (solid line). 
Down: Exponent a of \ oc T a as a function of InT for (c) 
MgTiOBOs and (d) MgMnB 2 5 . 



ture dependence of the exponent a(T) for both sam- 
ples. Furthermore, we fit the experimental data of Fig 
Ilia), using l/x(T) = Tln 2 (Q. /T) (solid line). Both 
Figs, ^a) and (c) indicate that the susceptibility co- 
incides exactly with the RSP model and a follows the 
same tendency previously predicted by numerical calcu- 
lations |2l|. In MgMnB205, previous assessments and 
the inset of Fig. IHb) suggest a power law behavior for 
x(T) with a constant a ~ 0.55. Within a more detailed 
analysis, shown in Fig.^d) we see an unequivocal slow 
increase of the exponent a, followed by a constant regime 
at intermediate temperatures. Although a is not con- 
stant in the whole temperature interval, as expected for a 
Griffiths phase (GP), its increase with T is clearly incon- 
sistent with a phase governed by an infinite randomness 
fixed point or RSP. However, for temperatures higher 
than 7 K, a is constant (up to 20 K), and this strongly 
supports the existence of a GP in this system. In fact 
the variation of a at low T may be related to a freezing 
of the Mn moments, as suggested by a low t emp erature 
anomaly in the specific heat of this material |ll| . 

Wc further investigate these two systems by comparing 
other independent measurements, such as magnetization 
and ac susceptibility as a function of a magnetic field, 
as shown in in Fig. [21 From Fig. Ufa) and Fig. |2Ib) we 
see that the MgTiOBOs data always present logarithmic 
corrections and the magnetization follows M oc ln(£>), 
as expected for a RSP |2(|. On the other hand, for 
MgMnB 2 5 both x(B) and M(B) follow a power law 
behavior with exponents a ~ 0.55 and 1 — a ~ 0.45 re- 
spectively. Such behavior is expected for systems in a GP. 
Finally, for the MgTiOBOs, we also analyze xa.c x T for 
different applied fields B (Fig.^c) ): the RSP is robust 
to applied fields bellow 3T even at temperatures up to 
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FIG. 2: (a) Magnetization normalized by the saturation mag- 
netization M a =g[iBNs and (b) a.c magnetic susceptibility as 
a function of the applied magnetic field for MgTiOBOa (open 
squares) and MgMnE^Os (open circles). The solid lines rep- 
resent a power law fitting for the MgMnl^Os data and a 
logarithmic fitting for the MgTiOBOa data, c) a.c magnetic 
susceptibility as a function of temperature for different values 
of an applied d.c field (MgTiOBOa). The solid line indicates 
a Curie-law fitting at high temperatures. 
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FIG. 3: (a) Experimental data of magnetic susceptibility for 
MgTiOBOa (closed circles) and the limit for the EW (solid 
line) . The system presents entanglement below the solid line. 
On panel (b) we show the same analysis for MgMnB2 0s. The 
insets show the quantity E as defined in eq. Q which quan- 
tifies the entanglement detected by the EW as a function of 
temperature. 



100K, where the susceptibility is not Curie-like. However, 
the RSP characteristics disappear at high temperatures 
once the applied field is around 3T with the appearance 
of a Curie-like behavior. 

Once established that MgTiOBOa is in a RSP, we can 
expect the system to be entangled. Theoretically, the 
entanglement can be estimated by calculating the Von 
Neumann entanglement entropy of a subsystem A of the 
spin chain, with respect to a subsystem B. This quantity 
can be defined as S = — Ttpa^Pa- For a subsystem 
with length x embedded in an infinite system, the entan- 
glement entropy for a Heisenberg chain in a RSP is given 
by S = ^g^ln(x), as previously calculated by means of 
a RSRG [10( and further confirmed by numerical calcula- 
tions |l5| . From an experimental point of view, it is not 
possible to use the entanglement entropy to quantify the 
entanglement. Entanglement witnesses (EW) have been 
proposed and applied as an attempt to detect and quan- 
tify entanglement experimentally. The main advantage 
of EW is that they are obscrvables and their mean value 
can be directly measured. 

The use of magnetic susceptibility as an EW is based 
on the violation of local uncertainty relations ■ The 
uncertainty of an operator Ai for a given quantum state 
is defined as AAf = (A 2 ) — (Ai) 2 , the statistical variance 
of the measurement outcomes . The uncertainty AA 2 can 
only be zero if the quantum state is an cigenstatc of Ai . A 
quantum state with zero uncertainty in all the properties 
Ai, must be a simultaneous eigenstate of these operators. 
If the simultaneous eigenstate does not exist, there must 
be a lower limit for the sum of the uncertainties. We can 

illustrate this concept for spins where s = (JV — l)/2: for 

-2 ~ 2 ~ 2 

any given spin state, we have [S x +S y + S Z )|V») = s(s+ 
1)|V>)- On the other hand, the expectation values of Si 
defines a vector with maximal length equal to ±s. Using 
both relations, we obtain the inequality AS 2 + AS 2 + 

AS 2 = (S x 2 + S v 2 + S z 2 ) - ((^) 2 + (S y ) 2 + (S z ) 2 ) > s 



We can apply this relation to obtain a limit for the 
correlation of separable states. Exemplifying for two 
spins 1 and 2: for product states, the uncertainty of 
Si(l) + Si(2) is equal to the sum of the local uncer- 
tainties. On the other hand, maximally entangled states 
can have a total uncertainty of zero in all directions of 
Si(l) + Si(2), maximally violating the uncertainty rela- 
tion. In this case, the quantity E = 1— A ( g '( 1 H g '( 2 )) 
measures the amount of entan glem ent verified by the 
violation of local uncertainties |22L l23l |. For a macro- 
scopic system, we can generalize this quantity by us- 
ing the magnetic susceptibility, which can be written as 
X, = ^Mi = ^((M 2 ) - ( Ml ) 2 ). 

Following ref[l|, the entanglement can be measured by 
the quantity 

E=l-k B T( X * + % +X > ). (1) 

In our case, N is the total number of spins per gram. 
First, we analyze a specific limit: if there is no dc 
applied field and the system is isotropic, E = 1 — 
3kBTxz/(g 2 t L B 2 Ns). Since the studied samples have 
vanishing magnetic anisotropy, only one component of 
the susceptibility is needed to detect and quantify entan- 
glement. For Xz < 9 3i^x S ' t ne system is entangled, and 
E quantifies the entanglement verified by the EW. 

In Fig. [3] wc show the experimental data for 
MgTiOB0 3 and MgMnB 2 5 : both systems present en- 
tanglement, although MgTiOBOa is entangled up to 
higher temperatures. The quantity E has a similar be- 
havior as a function of temperature for both compounds, 
with a linear dependence on T for high temperatures. 

For an applied dc field in the z direction, a pair of spins 
1/2, where J = Si + S 2 , form a singlet (H = ISi ■ S 2 ) 
at low fields. As [H, BJ Z ] = 0, the field does not modify 
the eigenstates, changing only their energies. At a given 
field, the energy of the singlet is no longer the lowest 
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energy B c , the singlet breaks and the spins align with 
the field B. However, for the whole range of fields, the 
ground state of the system is such that spin variance is 
minimal: AJ x AJ y = \{J Z ). In this case, as the system is 
isotropic in the xy plane, we have A 2 J y = A 2 J x = \{J Z ) 
so A 2 J y + A 2 J X = (J z ). This approximation is valid for 
g^ssB <C ksT, which assures that other states, which 
do not have this property, are not populated. Similarly, 
the same approach holds for two pairs of spin 5/2 as 
shown in Figs. E{a) andEJb). As an illustration, we also 
consider a distribution of singlets, where the probability 
for interaction strength / follows a power law. As can be 
seen in panel (c) of the same figure, the approximation 
works well for high values of the magnetic field compared 
with the temperatures. Since both systems are in a phase 
where the spins form dimers (MgTiOBOa is in a RSP and 
MgMnB205 presents Griffiths singularities in a random 
dimer phase) we can use this approximation to study 
quantum chains. It is possible to re- write the EW as 



E = 1 



9VbNs 



Xz 



(2) 



which is valid only at high fields. We perform the neces- 
sary measurements and using Eq. ^ for B = and Eq. 
for high magnetic fields (gfiBsB > 6/cbT for MgMnF^Os 
and guBsB > 2k B T for MgTiOB0 3 ) we quantify the en- 
tanglement for both systems. In Fig. 0] we unambigu- 
ously show that the magnetic field can increase entangle- 
ment in quantum spin chains, as theoretically suggested 
vcdrall, saguial. In the insets, we extrapolate the behav- 
ior of E x B for higher field values; measurements were 
performed with fields up to 9 T. From this extrapolation, 
we see that even if a field of 9T is not high enough for 
the approximation made in eq0 the extrapolation shows 
that the suggested increase of entanglement for low fields 
is still valid although the amount could bc slightly over- 
estimated. 

In conclusion, by means of macroscopic magnetic mea- 
surements we fully characterize a random singlet phase 
in a low dimensional spin system and for the first time, 
it was possible to distinguish this phase from a Griffiths 
phase. We use a novel analysis where the magnetic sus- 
ceptibility plays the role of an entanglement witness and 
measure the entanglement in two different spin systems 
as a function of temperature and magnetic field. We be- 
lieve that both types of analysis presented here can be 
used to experimentaly characterize the phase diagram of 
low dimensional systems. 
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FIG. 4: Up: Expectation value of the z component of the total 
spin {J z ), the sum of the total spin variance in three directions 
EiA 2 Ji, the spin variance in the z direction A 2 J Z and the 
sum A 2 J Z + {J z ) for (a) a pair of spins S=l/2, (b) a pair 
of spin S=5/2 and (c) pairs of spin S=l/2 interaction with 
a random interaction U (power-law distribution and T=0.05) 
as a function of magnetic field. This is normalized by the 
exchange interaction and in the random case by the cutoff of 
the distribution.. Down: Experimental data for E using eq. [5] 
for B for MgTiOB0 3 and MgMnB 2 5 . The insets show 
the extrapolation of E for very high values of B. 
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